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Abstract By analyzing the Chandra ACIS S3 data we studied 
the X-ray properties of the low-mass and high-mass X-ray binary 
populations in the nearby spiral galaxy NGC 5055. A total of 43 
X-ray point sources were detected within the 2 effective radii, with 
31 sources located on the disk and the rest 12 sources in the bulge. 
The resolved point sources dominate the total X-ray emission of 
the galaxy by accounting for about 80% of the total counts in 0.3- 
10 keV. By carrying out the spectral fittings we calculated the 0.3- 
10.0 keV luminosities of all the detected X-ray point sources and 
found that they span a wide range from a few 10 37 erg s _1 to over 
10 39 erg s -1 . After compensating for the incompleteness at the 
low luminosity end, we find that the corrected XLF of the bulge 
population is well fitted with a broken power-law model with a 

break at 1.57iolio x 1038 er § s -1 > while the P r °f>l e of the disk 
population's XLF agrees with a single power-law distribution with 
a slope of 0.931q Qg. The disk population is significantly richer at 
>2 x 10 38 erg s _1 than the bulge population, inferring that the 
disk may have undergone relatively recent, strong starbursts that 
significantly increased the HMXB population, although ongoing 
starbursts are also observed in the nuclear region. Similar XLF 
profiles of the bulge and disk populations were found in M81. 
However, in most other spiral galaxies different patterns of the 
spatial variation of the XLF profiles from the bulge to the disk 
have been observed, indicating that the star formation and evolu- 
tion history may be more complex than we have expected. 
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1 INTRODUCTION 

The study of the X-ray properties of the bright point sources in spiral galaxies, 
most of which are low-mass X-ray binaries (LMXBs) and high-mass X-ray bi- 
naries (HMXBs), may provide us with valuable observational constraints on the 
star formation and evolution history in the disk and bulge of the host galaxy. With 
the superb high spatial resolution of the Chandra X-Ray Observatory, a large 
population of X-ray point sources have been resolved for the first time in some 
nearby spiral galaxies. In a few specific cases, distinct X-ray characteristics have 
been revealed for the sources located in the bulge and those in the disk, which 
can be interpreted as the evidence for the spatial variation of stellar population 
composition that reflects the differences in star formation history (e.g., Tennant 
et al. 2001; Kong et al. 2002; Soria & Kong 2002). In this paper we present a 
Chandra study of the point sources in the nearby starburst galaxy NGC 5055 
(M63; SA(rs)bc), which is one of the prototype Arm Class 3 flocculent galaxies 
(Elmegreen & Elmegreen 1987) that shows regular, two-arm spiral structure to 
a radius of 4.0 kpc in the near-infrared band. The optical center of the galaxy is 
identified at RA=13hl5m49.25s Dec=+42d01m49.3s (J2000; Maoz et al. 1996). 
The LINER nucleus is U V bright and is surrounded by luminous young star clus- 
ters, showing clear stellar absorption signatures (Maoz et al. 1998; Leitherer et 
al. 2002). The inclination and position angles of the galaxy are deduced to be 
58° and 103°, respectively (Garcia-Gomez & Athanassoula 1991). We organize 
the paper as follows. In §2, we describe the observation and data reduction. In 
§3, we present the imaging analysis. In §4, we investigate the X-ray properties 
of the detected point sources, which includes the temporal and spectral analy- 
sis as well as the calculations of hardness ratios and X-ray luminosity functions 
(XLFs). Finally, we discuss and summarize the results in §5 and §6, respectively. 
Throughout the paper, we quote errors at the 90% confidence level unless men- 
tioned otherwise. We adopt a distance of 7.2 Mpc to NGC 5055 (Michele & 
Serra 1997), which is consistent with the distance calculated from its redshift 
(z = 0.001681) if cosmological parameters Hq = 70 km s _1 Mpc -1 , il m = 0.3 
and £1a = 0.7 are adopted. At this distance 1' corresponds to about 2.1 kpc. 

2 OBSERVATION AND DATA REDUCTION 

NGC 5055 was observed with Chandra with the CCD 0, 1, 2, 3, 6 and 7 of 
the Chandra Advanced CCD Imaging Spectrometer (ACIS) in operation on two 
separate occasions, which started on August 27, 2001 for a total exposure of 28.4 
ks, and on April 15, 2001 for a short duration of 2.4 ks, respectively. The events 
were telemetered in faint mode and the data were collected with a frame time of 
3.2 s. The CCD temperature was set at —120 °C. In the August 27 observation 
from which the data of this work was drawn, the center of galaxy was positioned 
on the ACIS S3 chip (CCD 7) with an offset of 31" from the nominal pointing for 
the S3 chip, so nearly all the emission of the galaxy was covered by the S3 chip. 
In the analysis that follows, we used the CIAO software version 2.3 to process 
the data extracted from the S3 chip only. In order to apply the latest calibration, 
we started with the Level- 1 data. We kept events with ASCA grades 0, 2, 3, 4 and 
6, and excluded bad pixels, bad columns, and columns adjacent to bad columns 
and node boundaries. In order to identify occasional periods of high background, 
we extracted the lightcurve of the source-free regions on the S3 chip in 2.5-7.0 
keV where the background flares are expected to be most visible. We found that 
up to about 5% of the total exposure time was affected by high background flares. 
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By excluding the contaminated intervals we obtained a clean exposure of 27.0 ks 
for the analysis. 

3 X-RAY IMAGE 

In Figure la we plot the Chandra S3 image of NGC 5055 in 0.3-10.0 keV in 
logarithmic scale. The image has been corrected for both exposure and back- 
ground, and has been smoothed by using a minimum signal-to-noise ratio of 

3 and a maximum signal-to-noise ratio of 5 per beam. We find that the X-ray 
emission from the galaxy is dominated by a large population of X-ray point 
sources, which contribute about 80% of the total counts of the galaxy in 0.3- 
10.0 keV. The diffuse X-ray emission is nearly symmetric within about 20", and 
is slightly elongated in east-west direction in outer regions. No diffuse X-rays is 
significantly detected beyond ~ 33", which is approximate the size of the bulge 
(~ 35"; Baggett et al. 1998). We find that the diffuse X-ray emission is peaked 
at RA=13hl5m49.3s DEC=+42d01m45.5s (J2000) where a bright nuclear X-ray 
point source is detected. The position of the X-ray peak also coincides with the 
optical and infrared centers of the galaxy to within 0.5". In Figure lb we show 
the optical image drawn from the Digital Sky Survey (DSS) in linear scale, on 
which the locations of the 43 detected X-ray sources are marked with circles for 
comparison (§4.1). 

4 X-RAY POINT SOURCES 

4.1 Detections 

We detected X-Ray point sources on the ACIS S3 image using the CIAO tool 
celldetect with a signal-to-noise threshold of 3. We first restricted the detections 
in 0.3-10.0 keV, and then crosschecked the results in 0.5-7.0 keV. We also have 
crosschecked the results by using the wavelet-based CIAO tool wavdetect and by 
visional-examination. On the whole S3 CCD 46 sources are detected, of which 
43 sources are located within the 2 effective radii (1 effective radius = \R e = 
1.6'; Thornley 1996). In this work we focus our study on these 43 sources only. 
The spatial distribution of the 43 sources show a clear concentration towards 
the galaxy center, inferring that most of them are physically associated with the 
galaxy. Based on the results of the deep Chandra observations of the blank fields 
(Mushotzky et al. 2000), we performed Monte-Carlo simulations and found that 
only few (<5) of the 43 point sources may be unrelated background sources. 
By overlaying the positions of these sources on the DSS image (Fig. lb), we 
find that 3 1 of the 43 sources are detected on the disk with a tendency to reside 
on the spiral arms. The rest 12 sources are detected in the bulge, of which up 
to about 1 source may be a disk population source that is misclassified into the 
bulge population due to the project effect. At the X-ray peak a nuclear source 
is detected, however it does not have the highest count rate among the detected 
sources. We list the properties of the 43 sources in Table 1, where we sort them 
in the order of increasing projected distance from the center of the galaxy. 

4.2 Temporal Variabilities 

After removing the intervals of strong background flares, we extracted the 0.3- 
10.0 keV lightcurves of the detected sources that each has more than 60 counts 
in total. We calculated the Kolmogoroff-Smirov (K-S) statistic for each of the 
lightcurves against the null hypothesis that the count rate of the source is tem- 
porally invariant over the effective exposure time; if the source is temporally 
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invariant, the cumulative fraction of the count is a diagonal from to 1 . We find 
that Src 22 (243 counts), Src 35 (600 counts) and Src 38 (102 counts) show 
significant evidence for temporal variations on the 90% confidence level (Fig. 
3). This is not likely to be caused by the variations of the local background, 
since the K-S test gives negative results in the background variability. For the 
central source (Src 1), the temporal variability is less significant in terms of the 
K-S test, possibly because K-S tests are most sensitive around the median value 
of the independent variable. However, by calculating the variability parameter 

S = (/max - /min)/^o^ max + er 2 -^ = 1.1, where / max and / min are the maxi- 
mum and minimum count rates, respectively, and c/ max and u/ min are the corre- 
sponding errors, we find that Src 1 is marginally variable on timescales of 1 ~ 5 
hr during which its count rate changed by about 50%. 

4.3 Hardness Ratios 

Since the background-corrected count rates of most of the resolved X-ray point 
sources are low, which makes it impossible to carry out spectral analysis for 
each source, we turned to study the hardness ratios of all the resolved sources 
that are defined as H21 = (M - S)/(M + S) and if 31 = (H — S)/(H + 
S), where S, M and if are the background-corrected counts extracted in 0.3- 
1.0 keV (S), 1.0-2.0 keV (M) and 2.0-10.0 keV (if), respectively. The same 
approach has been adopted in earlier works (e.g., Sarazin et al. 2000). We list the 
calculated hardness ratios and the lcr errors in column 7 and 8 of Table 1, and plot 
if 31 against if 21 in Figure 4. In the figure we show the predicted hardness ratio 
distributions for an absorbed power-law model by adopting column densities of 
1.31 x 10 20 cm" 2 (the Galactic value; Dickey & Lockman 1990) and 3.93 x 
10 20 cm~ 2 , and photon indices of T = 0.0, 1.0, 2.0, 3.0 and 4.0, as well as 
the hardness ratios for an absorbed blackbody model by adopting the Galactic 
absorption and temperatures of kT = 0.5, 0.4, 0.3, 0.2 and 0.1 keV. 

We find that the distribution of the colors is similar to that of M3 1 (Kong et 
al. 2002) and the nearly face-on, gas-rich spiral galaxy M83 (Soria & Wu 2003). 
Most of the sources lie in a broad diagonal band extending from (if 21, if 31) = 
(—1, —1) to (1,1). One source (Src 19) that is located at a moderate distance 
(1.1') to the galaxy center has hardness ratios of about (— 1, —1) and is thus 
identified as a supersoft source (SSS). If the selection criterions of Supper et al. 
(1997) and Kahabka (1999) are applied, 5 other sources (Src 6, 7, 9, 13 and 14) 
can be identified as candidates of SSSs as well since they satisfy ff31 + uh31 < 
— 1 and ff21 < 0, or if21 + 07/21 < —0.8. If we slacken the selection limit to 
ff21 < 0.5 and ff31 < 0.5 (e.g. Swartz et al. 2002), Src 12 and 15 can also 
be classified as SSS candidates. As is shown in Figure 4 the hardness ratios of 
these supersoft sources can approximately be described with an absorbed black- 
body model with the temperatures ranging from about 0.2 to 0.3 keV. Src 31 
has a hardness ratio of (1,1), which appears to have been heavily absorbed by 
a column density significantly larger than the Galactic value. Since Src 31 is 
located far away from the center of the galaxy (d > 110") and its hardness ratios 
agree with those expected by an absorbed power-law model with an absorption 
of > 0.4 x 10 22 cm~ 2 and a photon index of T ~ 0, it is probably an unrelated 
background AGN. 

For comparison, within the 2f? c region the measured mean hardness ratios 
for the total X-ray emission of the galaxy are (if 21, ff31) = (-0.28, 0.18). For 
the same region, the mean hardness ratios are (if 21, ff31) = (-0.14, 0.13) for 
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all the resolved X-ray point sources and (—0.46, 0.23) for the unresolved diffuse 
emission. 

4.4 Spectral Analysis 

We first extracted and studied the cumulative spectrum of all the 42 off-center 
point sources resolved within 2R C . The backgrounds were extracted from a care- 
fully selected annulus around each source. By applying the latest CALDB we 
have corrected for the charge transfer inefficiency (CTI) and the continuous 
degradation in the ACIS quantum efficiency, which is especially severe at lower 
energies. To avoid the effects of calibration uncertainties at lower energies and 
instrumental background at higher energies we restricted the spectral analysis to 
the 0.7-7.0 keV energy band. We found that the cumulative spectrum cannot be 
fitted with a single absorbed power-law model, if the absorption is fixed at the 
Galactic value. When the absorption was allowed to vary, however, the absorbed 
power-law model can give an acceptable fit (x 2 /dof = 152.7/116). The result- 
ing absorption and photon index are 0.14 ± 0.01 x 10 22 cm~ 2 and 2.05 ± 0.06, 
respectively. With these parameters the total flux of all the 42 off-center sources 
in 0.3-10.0 keV is calculated to be 7.79 x 10~ 13 erg s _1 cm~ 2 . We also have 
divided the resolved off-center sources into the disk and the bulge populations 
and studied their cumulative spectra in 0.7-7.0 keV. The cumulative spectrum of 
the disk population can be marginally fitted with an absorbed power-law model 
with an absorption of 0.27± 0.03 x 10 22 cm~ 2 and a photon index of 2.45 ± 0.08 
(X 2 /dof = 138.9/112). The spectrum of bulge population, on the other hand, 
can be well fitted with a model that consists of a power-law component and a 
black body component, both subjected to a common absorption that is consistent 
with the Galactic value (x 2 /dof = 60.1/55). The obtained photon index and 
temperature are T = 1.36 ± 0.13 and kT = 0.12 ± 0.01 keV, respectively. 

By excluding all the detected point sources we examined the spectrum of 
the diffuse emissions extracted in < 2R C that includes the contributions from 
both the unresolved point sources and the inter-stellar medium (ISM). The back- 
ground was extracted in a source-free boundary region on the S3 chip as far 
away as possible from the galaxy. We first attempted to fit the spectrum with an 
absorbed power-law or an absorbed MEKAL model, but neither of them gives 
an acceptable fit to the data, unless the absorption is allowed to increase to un- 
reasonable, physically meaningless values. Thus we attempted to apply a model 
that consists of a MEKAL component to represent the emission of the hot ISM, 
and a power-law component to represent the contribution of X-ray binaries, with 
both components subjected to a common absorption. Since the abundance of 
the MEKAL component is poorly constrained, we tentatively fixed it to 0.1 Zq, 
which is obtained by Tyler et al. (2004) for the central region of NGC 5055. 
We found that when the absorption is fixed to the Galactic value, the model pro- 
vides us with an acceptable fit (\ 2 /dof = 63.4/51), with a gas temperature of 
0.30 ± 0.02 keV and a photon index of 2.47 ± 0.26 . The total flux of the emis- 
sion in 0.3-10.0 keV is 7.24 x 10~ 13 erg s _1 cm~ 2 , of which about 54% can be 
ascribed to the power-law component. 

Of the 43 X-ray sources resolved within 2R C , 10 sources (Src 1, 5, 8, 9, 
16, 18, 22, 25, 27 and 35) have a total counts larger than 100. We extracted the 
individual spectra of these sources and the corresponding background spectra 
in annular regions adjacent to where the source spectra were extracted. We fitted 
each spectrum with an absorbed power-law model and/or an absorbed multicolor 
disk blackbody model (DBB), and list the results in Table 2. The spectrum of the 
center source (Src 1) is well fitted by an absorbed power-law model with a photon 
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index of 1.61 ±0.15, which is typical for those of supermassive black holes in ac- 
tive galaxies. The 0.3-10.0 keV luminosity of Src 1 corrected for the absorption 
is calculated to be 3.16xl0 38 erg s" 1 , which is lower than that of most low- 
luminosity AGNs. The spectrum of the brightest source Src 35, which showed 
significant temporal variations during the observation (§4.2), can be fitted with an 
absorbed DBB model with an inner disk temperature of kT in = 0.49 ± 0.04 keV. 
This source is a candidate ultra-luminous X-ray source (ULX) since its 0.3-10.0 
keV luminosity is 1.25 x 10 39 erg s _1 . 

According to the previous studies of the SSSs in spiral galaxies (e.g., M101, 
Pence et al. 2001; M31, Kong et al. 2002; M83, Soria & Wu 2003), SSSs are 
probably white dwarfs fueled by accretion from their low-mass companions. 
Indeed we find that the spectrum of the SSS candidate Src 9 can be fitted with 
an absorbed DBB model with the inner temperature of 0.15 ± 0.01 keV, which 
is consistent with that of an accreting white dwarf. 

4.5 X-Ray Luminosity Functions of the Resolved Off-Center Point Sources 

All Off-Center Sources Detected within 2R C 

Assuming that all the off-center X-ray point sources resolved within 2R C 
are located at the distance of NGC 5055, we calculated their X-ray luminosi- 
ties in 0.3-10.0 keV using the best-fit spectral parameters for their cumulative 
spectrum. The conversion factor of the counts is 2.13 x 10 36 erg cts -1 , and 
the resulting luminosities range from 2.1 x 10 37 to 1.25 x 10 39 erg s _1 . With 
these we construct the XLF and illustrate it in Figure 5a. Since the detection of 
the point sources is not complete at the faint end of the luminosity function, by 
adopting a method similar to that utilized in, e.g., Kim and Fabbiano (2004) and 
Xu et al. (2005) we ran Monte-Carlo simulations to create fake point sources on 
the S3 image of NGC 5055 in 0.3-10.0 keV. In the simulations we created fake 
sources by using the MARX package (Wise et al. 1997), and assumed that the 
radial distribution of the fake sources at any given luminosity follows the r 1 / 4 
law (de Vaucouleurs 1948). At a given luminosity, we determined how many of 
the fake sources can be detected with the same technique used in §4.1. In such a 
way we corrected both the observed XLF and the background for the unresolved 
sources. The corrected XLF are also shown in Figure 5b. 

We fitted both the uncorrected and corrected cumulative XLFs with the soft- 
ware Sherpa by using either a single power-law profile or a broken power-law 
profile 



where L/L 38 is the 0.3-10.0 keV luminosity in units of 10 38 erg s _1 , and a; and 
ah are the slope indices for the lower and higher luminosity ends, respectively. 
We find that for the uncorrected XLF, the fittings with the single power-law 
model can be immediately rejected (x 2 /dof = 150.2/21). It overestimates the 
data at the high luminosities and underestimates the data at the low luminosities. 
The use of the broken power-law model, on the other hand, can significantly 
improve the fittings and provide us with an acceptable fit (x 2 /dof = 15.8/19). 
The best-fit parameters at the 90% confidence level are L b = 2. 53^45 x 10 38 
erg s"\ a t = 0.61±o;°,3 and a h = 1.96:^8. For the XLF corrected for the 
effect of incompleteness at the faint end of the XLF, the single power-law model 
also cannot give an acceptable fit (x 2 /dof — 125.8/22), while the broken 
power-law model can improve the fittings significantly (x 2 /dof — 23.1/20) 
with the best-fit parameters L b = 3.11±g;|| x 10 38 erg s"\ a L = 0.71 ± 0.05 
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and a h = 2.36tg;g. 

Bulge and Disk Populations 

We studied the XLFs of the bulge and disk population sources (Fig. 5) 
and found that the single power-law model is inadequate to describe both the 
uncorrected XLFs of the bulge population (x 2 /dof = 68.9/9) and the disk 
population (x 2 /dof = 30.3/19). The broken power-law model, however, gives 
a good fit to the data with L b = 1.68^.22 x 10 38 erg s" 1 , a t = 0.31 ± 0.12 
and a h = 3.95±\H for the bulge population (x 2 /dof = 7.4/7), and 
L b = 3.1lt^8 x 10 38 erg s"\ a t = 0.71 ± 0.05 and a h = 2.36±g;t| for 
the disk population (x 2 /dof = 23.1/20). We then corrected the XLFs for the 
incompleteness at the low energies using the method as is described above for 
all the resolved off-center sources. For the bulge population, a broken power-law 
is still needed to describe the XLF with L b = 1.57±o2o x 10 38 erg s~\ 
a t = 0.38±g;|g and a h = 2.28t*; 3 ° (x 2 /dof = 7 A/7). The corrected XLF 
of the disk population is found to be nicely consistent with a single power-law 
model with a slope of a = 0.93 tom (X 2 /dof = 16.3/19). 

We have crosschecked our results on the XLF profiles of all the resolved off- 
center sources and disk population by excluding Src 35, the brightest source that 
may bias the fittings. We found that within the errors the best-fit parameters are 
consistent with those obtained with Src 35 included. 

5 DISCUSSION 

We calculated the 0.3-10.0 keV luminosities of 43 X-ray point sources detected 
within 2R C of NGC 5055 and found that they span a wide range from about 
2.1 x 10 37 erg s _1 to 1.25 x 10 39 erg s _1 , which is typical for spiral galaxies 
(M81; Tennant et al. 2001, NGC 1637; Immler et al. 2003) and early-type galax- 
ies (Xu et al. 2005 and references therein) at similar distances. After compensat- 
ing for the incompleteness at the low luminosity end, we find that the corrected 
XLF of the bulge population is well fitted with a broken power-law model, while 
the profile of the disk population's XLF satisfies a single power-law distribution. 
The disk population is significantly richer at the high-luminosity end (>2 x 10 38 
erg s _1 ) than the bulge population, inferring that the star formation history of 
the bulge is distinct from that of the disk. In other words, the disk may have 
undergone recent, strong starbursts that significantly increased the HMXB pop- 
ulation, although ongoing starbursts are also observed in the nuclear region of 
NGC 5055. Quite similar phenomenon has been found in M8 1 (also identified as 
a LINER or Syl.8) by Tennant et al. (2001), who reported that the XLF of the 
bulge population exhibits a break at ~ 4 x 10 37 erg s _1 , and the XLF of the disk 
population can be fitted with a single power-law. 

Different patterns of the spatial variation of the XLF profiles from the bulge 
to the disk are observed in other spiral galaxies. In M83, a SAB(s)c galaxy that 
hosts a starburst nucleus (Soira & Wu 2003), the profile of the bulge XLF is con- 
sistent with a single power-law with a slope of 0.8, while a break appears on the 
XLF of the disk population at about 8 x 10 37 erg s _1 . This might suggest that, un- 
like in NGC 5055 and M8 1 , there is no strong starbursts in M83 on the disk in the 
recent past, and the HMXBs dominate the source population only in the nuclear 
starburst regions, rather than over the entire disk. In M31 Primini et al. (1993), 
Shirey et al. (2001), Kaaret (2002) and Kong et al. (2002) found that the cumula- 
tive XLF of all the detected point sources has a distinct break at 2 — 3 x 10 37 erg 
s _1 , a slope of 1.36 ±0.33 for the high luminosity end and a slope of 0.47 ±0.09 
for the low luminosity end. Kong et al. (2002) also studied the cumulative XLFs 
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of the point sources in the inner bulge, outer bulge and disk regions. They found 
that there is a spatial variation in the XLF profiles; both the break luminosity and 
the slope of the cumulative XLF increase outwards monotonously from the inner 
bulge, where the LMXBs dominate the X-ray emission of point sources, to the 
disk. In a few other spiral galaxies where the XLFs have been studied in detail, 
however, the report of a break on the XLF is not available in literatures. These 
includes M101 (Pence et al. 2001), IC342 (Bauer et al. 2003; Kong 2003), NGC 
891 (Temple et al. 2005), NGC 1637 (Immler et al. 2003) and M51 (Terashima 
& Wilson 2004), for which the slopes of the XLFs range from about 0.5 to about 
1. 

The diversity of the XLF profiles of the X-ray sources in spiral galaxies in- 
dicates that the star formation and evolution history may be more complex than 
we have expected, so that further careful multi-band investigations are needed. In 
addition, we suggest that, since typically only 50- 150 X-ray point sources are de- 
tected per galaxy, the uncertainty introduced by the small number statistics may 
bias our conclusions on the XLF profile to a certain extent. In Xu et al. (2005), 
by performing Monte-Carlo simulations we showed that even if there is an uni- 
versal break on the XLFs of early-type galaxies, the statistical errors definitely 
preclude us from measuring it correctly. The same result can be applied to the 
measurements of the XLF profiles of spiral galaxies. 

6 SUMMARY 

We detected a total of 43 X-ray point sources (12 in the bulge and 31 on the disk) 
within the 2 effective radii of NGC 5055, whose 0.3-10.0 keV luminosities range 
from a few 10 37 erg s _1 to over 10 39 erg s _1 . After compensating for the incom- 
pleteness at the low luminosity end, the corrected XLF of the bulge population 
is well fitted with a broken power-law model (L b = 1-57±q%1 x 10 38 erg s _1 ), 
while the XLF of the disk population satisfies a single power-law distribution 
(a = 0.93^Q'Qg). The disk population is significantly richer at >2 x 10 38 erg s _1 
than the bulge population, inferring that the disk may have undergone relatively 
recent, strong starbursts that significantly increased the HMXB population, al- 
though ongoing starbursts are also observed in the nuclear region. This is similar 
to what has been found in M8 1 . In most other spiral galaxies, however, different 
patterns of the spatial variation of the XLF profiles from the bulge to the disk are 
observed, indicating that the star formation and evolution history may be more 
complex than we have expected. 
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Fig. 1 (a): X-ray image of the central 3.6' x 3.6' region of NGC 5055 
in 0.3-10 keV in logarithmic scale, which has been smoothed with a 
minimum significance of 3 and a maximum significance of 5, and has 
been corrected for both exposure and background, (b): DSS blue image 
with the field of view of (a) and the bulge marked with a box and a 
circle, respectively. On both images the detected X-ray point sources are 
marked with small circles. 
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Fig. 2 Background-corrected lightcurves of Src 22, Src 35 and Src 38. 
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Fig. 3 Results of Kolmogoroff-Smirov test for the lightcurves of Src 22, 
Src 35 and Src 38. 
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a The columns are arranged as follows: (1) source number; (2)-(3): right ascension and declination 
(J2000); (4) count rate and its error; (5) projected distance d to the center of the galaxy; (6) intrinsic 
X-ray luminosity L x , assuming that the source is located at the distance of NGC 5055 and only 
subjected to the Galactic absorption (1.31 X 10 20 cm" 2 ; Dickey & Lockman 1990); (7)-(8): 
hardness ratios and the 1<t errors (§4.3); and (9) notes. 

b - c H 2 i = (M-S)/(M+S), H :n = (H-S)/(H+S) 

d Abbreviations are b: source in bulge; c: central source; s: candidate of supersoft source; u: 
ultra luminous x-ray source; and v: source shows temporal variability. 
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Table 2 Spectral Fittings of the Brightest Sources a 
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1 Errors are quoted at the 90% confidence level. 
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Fig. 4 Color-color diagram of the resolved point sources. Super soft 
sources (SSSs) and the hardest source are marked with open circles and 
an open star, respectively. Model predictions are shown with dashed and 
dotted lines. From Top of the figure downwards, open rectangles on the 
lines indicate the values of photon indices of T = 0.0, 1.0, 2.0, 3.0 and 
4.0 for power-law models, and temperature of kT = 0.5, 0.4, 0.3, 0.2 
and 0. 1 keV for blackbody model, respectively. 



14 



Luo et al. 



A 
55 



A 
S5 



-i 1 — i — i — i i i 1 1 

integrated 



-i — i — i — i i i i | 




(a) Uncorrected 

\ 1 i i i i 1 1 1 1 h 



integrated 



(b) Corrected 



i i i i i 1 1 




-1 I I I I I I I I I L 



10 100 
L x [0.3-10 keV] (10 37 erg s -1 ) 



Fig. 5 XLFs of the detected off-center X-ray point sources uncorrected 
(a) and corrected (b) for the incompleteness at the low luminosity end. 
Models are plotted as dashed lines. A single power-law model is suffi- 
cient to describe the corrected XLF of the disk population, while in other 
cases a broken power-law model is needed (§4.5). 



